Exact results in a slave boson saddle point approach 
for a strongly correlated electron model 
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We revisit the Kotliar-Ruckenstein (KR) slave boson saddle point evaluation for a two-site cor- 
related electron model. As the model can be solved analytically, it is possible to compare the KR 
saddle point results to the exact many particle levels. The considered two site cluster mimics an 
mfinite-(7 single-impurity Anderson model with a nearest neighbor Coulomb interaction: one site is 
strongly correlated with an infinite local Coulomb repulsion which hybridizes with the second site, 
on which the local Coulomb repulsion vanishes. Making use of the flexibility of the representation we 
introduce appropriate weight factors in the KR saddle point scheme. Ground state and all excitation 
levels agree with the exact diagonalization results. Thermodynamics and correlation functions may 
be recovered in a suitably renormalized saddle point evaluation. 
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The unusual properties of the high temperature super- 
conductors or the colossal magnetoresistance and orbital 
ordering in transition metal oxides^ are calling for new 
techniques and concepts to describe these phenomena. A 
formalism, devised to deal specifically with strongly cor- 
related electron systems, is provided by the slave boson 
approach pioneered by Barnes 2 - for the single impurity 
Anderson model (SIAM), and further extended to the 
Hubbard model by Kotliar and Ruckenstein.— The ap- 
proach implements a local decomposition of electronic 
excitations into charge and spin components which is 
achieved by the introduction of composite operators for 
all local (on-site) excitations. The composite operators 
separate into canonical bosonic and fermionic operators 
where, in the Barnes approach, the charge (spin) degrees 
of freedom are represented by bosonic (fermionic) op- 
erators. However the latter are enslaved in the sense 
that their respective number operators have to fulfill 
a local constraint. The original idea was to decouple 
spin and charge degrees of freedom. In the Kotliar- 
Ruckenstein (KR) scheme a bosonic mode is attributed 
to each type of excitation which allows to study the cor- 
related system in a saddle point approximation for all 
degrees of freedom.— 4£ This latter approach has been im- 
pressively successful when compared to numerical simula- 
tions: ground state energies^ and charge structure factors 
show excellent agreement f 7 - 

The slave boson approach has several intriguing prop- 
erties. Firstly, the approach is exact in the large de- 
generacy limit4>& Moreover, the paramagnetic mean- 
field solution reproduces the Gutzwiller approximation 
in the KR representation.— It obeys a variational prin- 
ciple in the limit of large spatial dimensions where the 
Gutzwiller approximation and the Gutzwiller wave func- 
tion are identical.— These formal properties signify that 
the approach captures characteristic features of strongly 
correlated electrons as the suppression of the quasi- 



particle weight and the Mott-Hubbard/Brinkman-Rice 
transition^ to an insulating state at half filling with in- 
creasing on-site Coulomb interaction. 

Secondly, as the fields are canonical within this ap- 
proach, long range correlations are more amenable to 
analyses. Moreover, one can easily introduce a long range 
Coulomb interaction. In a radial gauge representation of 
the slave bosons, long range Coulomb interactions can 
even be cast into a bilinear form.— In other approaches, 
which satisfactorily implement local correlations, long 
range Coulomb interactions are typically difficult to han- 
dle. 

The flexibility of the slave boson approach has con- 
tributed to its success in many fields from Kondo 
physic a 12 ' 13 and Kondo/ Anderson lattice d 14 ' 15 ' 16 to inter- 
faces of correlated electronic systems . 17 ' 18 It separates in 
a straightforward way the high and low energy scales. By 
now the respective mean-field evaluations are well doc- 
umented (see, e.g., Refs. |f>. 19.20.21.22]), and a number 
of fluctuation calculations have been performed,— i 23 i 24 i 25 
even though the choice of the proper framework has been 
intensely debated) 4 - ' 11 ' 23 ' 26 ' 27 In the case of the single im- 
purity Anderson model (SIAM) one rather resorts to di- 
agrammatic approaches which include slave boson tech- 
niques (for a recent reference see Ref. [281]). 

It was doubted that the low energy physics is imple- 
mented consistently in the mean field evaluations because 
the decomposition is a local scheme. In the present work 
we want to access all energy scales for a specific model 
and compare the saddle point evaluation of the KR slave 
boson approach to exact results. Comparing exact and 
saddle point calculations analytically represents a diffi- 
cult task when handling a large system. Therefore, as a 
first step, we restrict our considerations to a two-site clus- 
ter which can be diagonalized exactly but still incorpo- 
rates characteristics of a strongly correlated electron sys- 
tem: the considered two site cluster presents a truncated 
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infinite-?/ single-impurity Anderson model with a nearest 
neighbor Coulomb interaction. One site is strongly corre- 
lated with an infinite local Coulomb repulsion which hy- 
bridizes with the second site, on which the local Coulomb 
repulsion vanishes. From the investigation of this model 
another benchmark of the slave boson approach is pro- 
vided here and, moreover, two essential questions may 
be addressed: how do thermodynamics and correlation 
functions compare between exact and saddle point eval- 
uation? How is the intersite Coulomb interaction appro- 
priately introduced in the KR slave boson approach? 

Interacting two-site cluster. We introduce the SIAM- 
type Hamiltonian 

H = £ ( e c4c CT + e d dtd a + V (cld a + h.c.)) 

+Ud^d^d^d^ + In&n c (1) 

which is defined on two sites. The operators cj. (c CT ) and 
d'l (d a ) represent the creation (annihilation) of "band" 
and "impurity" electrons, with spin projection a. The en- 
ergies e c and are the band and impurity energy levels, 
respectively, while the hybridization energy is V. Here U 
is the on-site repulsion, which is the largest energy scale 
in the model and it will be set to infinity as in the stan- 
dard infinite-?/ SIAM. Finally, Hi = Iridn c represents 
the non-local Coulomb interaction. 

This extended two-site SIAM may be solved exactly, 
either by diagonalization of the Hamiltonian, Eq. (TJ, or, 
in the Lagrangian language , 29 ' 30 through the exact evalu- 
ation of the path integrals representing the desired quan- 
tities within a slave boson representation in the radial 
gauged The two-particle basis for the Hamiltonian ma- 
trix consists of two singlet states and three triplet states. 
Indeed one finds a threefold degenerate eigenvalue of H 

= (e c + e d + I) (2) 

corresponding to the triplet states, and two non- 
degenerate eigenvalues controlled through A = e c — 

E^= ~ (3e c + e d + 1 ± V(A - If + 8V 2 ) (3) 

corresponding to the singlet states. 

KR slave boson scheme. The goal is now to formulate 
an effective theory which correctly reproduces the expec- 
tation value of the charge density and the charge and 
spin density correlation functions. The requirement to 
work with canonical fermionic or bosonic operators (or 
rather 'fields' in the Lagrangian language) leads to the 
slave boson scheme with which we intend to evaluate the 
free energy in the KR SBMF (slave boson mean field or, 
equivalently, saddle point) approximation. The first step 
in the scheme is to enlarge the Fock space for the impu- 
rity site by introducing auxiliary fermions f& and slave 
bosons and p&\ where e^, acts on an empty and 
Pa on a singly occupied impurity site. A field for double 



occupancy of the impurity site is omitted as we restrict 
the evaluation to infinite U. The electron operator on 
the impurity is then represented as a composite operator 

4 = Me, d„=JpJ a , (4) 

These fields are subject to three constraints: 

eh + ^2plp a = 1 

pIp, = flU ° =T, I (5) 

which are enforced by three Lagrange multipliers, de- 
noted by a and A CT , respectively. If strictly enforced, 
only empty and singly occupied states in the physical 
space remain. 

Now, in the second step of the setup, the KR technique 
renormalizes the coupling parameters of the Hamiltonian 
or Lagrangian. It thereby makes use of the freedom to 
choose different (operator) representations of the respec- 
tive coupling terms, which are all equivalent in the phys- 
ical subspace. This freedom originates from the decom- 
position of the electron field into auxiliary fields. Kotliar 
and Ruckenstein-. introduced z-factors into the hybridiza- 
tion term which then reads for the SIAM: 

H v = F^(4d CT +h.c.) 

= V£( c tz ff / ff +h.c.) . (6) 

a 

A straightforward choice would be z a = e' p a which di- 
rectly translates the relation Eq. (@]) within Hv- How- 
ever, the well established substitution is rather 

«a = e f (1 - pip.)'* (1 - eh - ptrp^)-* p a (7) 

which reproduces the U = spinless case and the 
Gutzwiller approximation^ in the saddle point evaluation 
and which led to the excellent agreement of the SBMF 
with numerical simulations in the case of the single band 
Hubbard model. We note that Hv has the same ma- 
trix elements in the physical subspace for any of these 
representations for z a . 

In the third step of the scheme the bosonic fields, in- 
cluding the Lagrange multipliers, are replaced by their 
respective saddle point values. Here we emphasize that 

is in fact the quasiparticle weight, as may be confirmed 
by straightforward evaluation of 9I q^ = 1 — -p-, where 

represents the impurity site self-energy. 
Before we discuss the KR SBMF results, the renormal- 
ization of the nonlocal Coulomb interaction / has to be 
introduced: 

Hi = In d n c = I y^(l - yjyo) n C;CT . (8) 

a 

The natural choice would be y a = e (spin independent), 
as (1 — eh) is the density operator. However, on the 
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saddle point level, we rather expect that y a = z\ should 
account for a more appropriate representation. As z\z a is 
related to the quasiparticle weight, \ — y\y a characterizes 
the incoherent part of an electronic excitation. It here 
refers to a local process, Eq. J8]). Correspondingly we 
define: 

Vt = e f (1 + £ - ylVoY^ (1 - e f e - ^- a V- a Y^ V a • (9) 

An infinitesimal convergence factor e > has been intro- 
duced to ensure the property y a — for e = 0, also on 
the saddle point level. This singular assignment of y a = 
may occur if the hybridization is suppressed either due 
to the formation of a triplet state or on account of the 
special choice V = 0. Otherwise e may be set to zero 
from the outset. Such a convergence factor is not neces- 
sary for z a in the hybridization, Eq. ([6]). We emphasize 
again that replacing In<in c by the expression on the rhs 
of Eq. ([8]) with y\ from Eq. J9]) does not hurt the cor- 
rectness of the representation, i.e., all matrix elements of 
Hi are unaffected in the physical subspace. Moreover, by 
extending the scheme developped in Refs. [29ll3oT | to the 
KR representation, it can be shown by direct evaluation 
of the path integral that the exact partition function is 
recovered when using Eqs. (0[9|). 

KR SBMF. The saddle point evaluation is now easily 
implemented: (i) the bosonic fields are replaced by real 
variables and (ii) these variables are found from the min- 
imization of the grand canonical potential. With (i) we 
identify the fermionic matrix as: 

"w- (iv Z)- m 

with E c , a — e c + (1 - y 2 ) I - \i and E {<(7 = e d + \ a - fi. 
Its eigenvalues are 

E±, a = \ {E^E Ua ±sj (E c , a - E iM f + 4z 2 V 2 ) (11) 

For convenience, we have introduced a joint chemical po- 
tential [i for all electrons (on the d- and the c-site) to 
control the filling. 

Part (ii) of the evaluation depends crucially on the 
number of electrons JV is the system. We first con- 
sider the two-electron case (JV — 2). The free energy F 
reads: F(N = 2,T) = -T£ CTr=± m(l + e X p(-AE r , CT ))- 
J2 a ^crPa~ a (l~ e2 ~J2aPa) + '2{i, where the first log-term 
is standard and the terms with the Lagrange multipliers 
A CT and a were generated through the constraints. The 
phases of the slave boson fields do not enter into the sad- 
dle point evaluation, so e and p a are now real numbers, 
as well as a and A CT . The quadratic term e 2 is to be iden- 
tified with the hole number x on the correlated site, i.e., 
x = e 2 . We take the temperature T in units of fcs and 
(3 = 1/T. For the limit T -> we have F(N = 2,0) = 

E CT E -,° ~ ^Pl + "(e 2 + J2a A ~ l ) + V Th e mini- 
mization of the ground state F(N = 2, 0) is cumbersome 



but straightforward. We obtain: 

F(N = 2 , T = 0) - 2e d = 

(A(l+x)+I(l-x)-2VVteVT^x) (12) 

Here, the hole number on the correlated site obeys the 
analytical relation (with v = — 1 and x_ = x): 

X v = - 1 + V — ; 13 

which displays a monotonous decay with increasing (A — 
I)/V. The ground state is paramagnetic and, with 
p 2 = p 2 _ a = p 2 , the constraint yields the correct rela- 
tion p 2 — (1 — x)/2 for the expectation value of single 
occupancy on the correlated site. Notably, these rela- 
tions, Eqs. (j 1 2 P 1 3(1 . are exact. This statement is con- 
sistent with the observation that the lowest two-particle 
state acquires the energy (E- t -\ + £_,]_) = E^ of the 
lower antiferromagnetic singlet, Eq. ([3]), for the saddle 
point solution of the bosonic fields. 

However, it is not only the lowest energy state in the 
two particle sector, which is solved exactly in the KR 
SBMF scheme. The approach also provides the exact 
solutions for the single particle (N = 1) and the three 
particle (N = 3) sectors, and moreover the triplet solu- 
tion for (JV = 2). These latter cases are not as surprising 
as the two-particle solution for the paramagnetic ground 
state. Nevertheless, these states will be important for the 
calculation of correlation functions. 

Here it suffices to emphasize that the KR approach 
was set up with the proviso that the spin polarized limit 
is correctly implemented through the weight factors z a . 3 
In fact, we find for JV = 1 the hol e expectat ion value 
on the correlated site x = — A/V A 2 + 4V 2 ) and the 
fermionic matrix Eq. (fTO)) reduces to the simple 2x2 ma- 
trix with entries e c and ed on the diagonal and V as the 
off-diagonal matrix elements. This represents correctly 
a two-level problem. The N — 2 triplet state and the 
JV = 3 case are exactly reproduced because these are local 
states with vanishing kinetic term Hv and the weight fac- 
tors y a for the nearest neighbor interaction I have been 
chosen suitably in Eq. ((9]). Truly, for the triplet state, 
we could verify that (J?-,t + E + ^) = E^> for E±^ from 
the corresponding saddle point solution. There, the solu- 
tion becomes the ground state if the model free energy is 
extended by a magnetic term, linear in the external field. 

The high energy singlet state presents, for the limit 
V — > 0, a local singlet with the two electrons on the un- 
correlated site. For V ^ 0, this singlet becomes non-local 
as the electronic states hybridize with the correlated site. 
It is important to observe that the KR SBMF scheme is 
flexible enough to accommodate also this second non- 
local two-particle state correctly. In order to gain the 
exact expression E + of Eq. (|3|) , the saddle point evalu- 
ation for the bosonic fields has to be retraced, however, 
with (-E?+,f + E + _i) for the two-particle state. We find 
for the hole number on the correlated site x — x + (see 



4 



Eq. (fT3t) and {E +A + E + . l ) = E^> . In other words, we 
here minimize the energy in the subspace of the Hilbert 
space that is orthogonal to the ground state. 

Spin Correlation function. The spin correlation func- 
tion involves transitions between states in the singlet and 
triplet spin sector which are separated by E^ — E_ , For 
the two-site model we can easily evaluate the spin corre- 
lation function ((f>± IS'+^S'-^)!^^) where <j>± refers to 
the singlet states and S_(0) = djc^ = S+(0) acts on the 
correlated (impurity) site. The singlet state ket-vector 
is: 

i4 s) ) = a (cfcpio) ± b (d{4 + 44)i°) ( 14 ) 

where the fermionic operators act on their vacuum state 
|0), and the normalization factors are a = y/x and b — 

v/|(l — x). In fact, the total spin operator for 
the two sites correctly yields the singlet state property 
£+°!?l0±^) = 0. The spin correlation function for the 

correlated site is, e.g., for T < E {t) - E^ and A - I > 
0, and for a chemical potential which fixes the particle 
number to two: 

(^le-^e-^dJe'^dl^W) 

The procedure for its evaluation in the KR technique 
is straightforward: the exponential Hamiltonian terms 
are taken in their respective mean field form which act 
on the states generated by the spin operators. As the 



saddle point evaluation delivers the exact energy levels, 
we identify: 

(S + (t)S.(0)) = l(l- X )^ Em - E ^ (15) 

where the partition function Z = er^ E - cancels the 
temperature factor. Similarly, the exact spin correlation 
function with a trace, which includes also the triplet and 
high energy singlet, is recovered for arbitrary tempera- 
ture values. 

It is the remarkable flexibility of the slave boson tech- 
nique which correctly renders certain limits of strongly 
correlated electron systems. This adaptability originates 
in the freedom to choose appropriate weight factors for 
the interaction terms in the projected field theory. For 
the two-site SIAM, the generic choice of this weight fac- 
tors already leads to the exact solution of the model. The 
choice is generic in the sense that it was introduced early 
on for the Hubbard model in the original work on the 
KR slave boson scheme.— It is striking to observe that 
even the inclusion of a non-local Coulomb repulsion in 
this model does not break the exactitude of the saddle 
point evaluation in this scheme. 
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